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We present the anisotropic optical conductivity of MgB2 between 0.1 and 3.7 eV at room temper- 
ature obtained on single crystals of different purity by the spectroscopic ellipsometry and reflectance 
measurements. The bare (unscreened) plasma frequency ui p is almost isotropic and equal to 6.3 eV, 
which contrasts some earlier reports of a very small value of ui p . The data suggests that the cr-bands 
are characterized by a stronger electron-phonon coupling Xtr but smaller impurity scattering 7i mp , 
compared to the 7r-bands. The optical response along the boron planes is marked by an intense 
interband transition at 2.6 eV, due to which the reflectivity plasma edges along the a- and c-axes 
are shifted with respect to each other. As a result, the sample spectacularly changes color from a 
blueish-silver to the yellow as the polarization is rotated from the in-plane direction towards the 
c-axis. The optical spectra are in good agreement with the published ab initio calculations. The 
remaining discrepancies can be explained by the relative shift of cr-bands and 7r-bands by about 0.2 
eV compared to the theoretical band structure, in agreement with the de Haas-van Alphen exper- 
iments. The widths of the Drude and the interband peaks are both very sensitive to the sample 
purity. 
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I. INTRODUCTION 

The discovery of superconductivity in MgBj 1 " caused 
an excitement in the physics community not only due 
to an unprecedentedly high T c (^40 K) for a 'conven- 
tional' electron-phonon superconductor, but also because 
it clearly shows the existence of two distinct gaps^; an in- 
triguing phenomenon that, although addressed theoret- 
ically, had been never observed before. Magnesium di- 
boride consists of graphite-like boron planes intercalated 
with Mg atoms. The metallic properties are determined 
by two distinct types of electronic bands: the strongly 
covalent almost 2D cr-bands formed by the hybridized 
sPxPy B orbitals and 3D 7r-bands made of p z orbitals. 
The holes in the cr-bands are strongly coupled to the 
in-plane bond-stretching phonon modes, giving rise to a 
high electron-phonon coupling constant. 3 A remarkable 
implication of such a strong conduction band disparity 
is a multigap superconductivity: a large gap value on 
cr-bands and a small one on 7r-bands. Notably, the com- 
munity has enjoyed a rapid advance in understanding 
the electronic structure and superconducting scenario of 
MgB2, in contrast to the lengthy siege of the high-T c 
problem in the cuprates. 

The early theoretical predictions of the electronic 
structure and the superconducting properties were 
soon confirmed by the isotope effect)^ angle-resolved 
photoemissionj& de Haas- van Alphen (dHvA)^ specific 
heati-iiLiiL tunneling measurements^ and, recently, by 
the inelastic X-ray scattering^ Even though the far- 
infrared experimentsi2ii^i£*i§iiLi£ have clearly shown 
the lowest of the superconducting gaps in agreement with 
the theory and other spectroscopic probes, the optical 



measurements so far demonstrated a rather poor repro- 
ducibility and equally poor consistency with the theo- 
retical electronic structure. The most controversial issue 
is value of the Drude plasma frequency lu p . While the 
theory predicts a high value of lu p of 7 eV, which cor- 
responds to about 1 conducting electron per unit cell, a 
number of g r0 upsi2ii2i22i2ii22i22i reported a much smaller- 
value of about 1.5 - 2.5 eV, corresponding to 0.15 elec- 
trons per cell. In Refl2fl in addition to the narrow Drude 
peak, a broad continuum, which could contain the miss- 
ing Drude spectral weight, was found below 1 eV. How- 
ever, the existence of such a continuum was not reliably 
verified by other groups. Another inconsistency is re- 
lated to the anisotropy of the plasma frequency. The 
calculations provide very close (within 5%) values of Lo p ^ a 
and ujp.c, which should be regarded, in fact, as a coinci- 
dence, because the plasma frequency of the carriers in the 
2D cr-bands is strongly anisotropic. A study of the op- 
tical anisotropy of magnesium diboride was undertaken 
by Fudamoto and Lee2£ by the comparison of the reflec- 
tivity spectra measured on a mosaic of ab-oriented crys- 
tals and on a polycrystalline sample. They observed the 
in-plane reflectivity plasma edge at about 16000 cm -1 . 
The authors also suggested that an additional structure 
in the reflection of a polycrystal is due to the c-axis 
plasma edge at 22000 cm" 1 . Thus the anisotropy ra- 
tio Wp. a /LUp^ c ~ 0.73 was deduced, which is in contrast 
with the theoretical prediction. 

The optically derived electron-phonon coupling con- 
stant Xtr was also a subject of debatei 19 ' 22 i 24 ' 25 Values of 
Xtr ranging from 0.13i£ to 1.5— were reported. It was 
pointed out in Ref. 25 that the determination of X tr re- 
lies strongly on the plasma frequency, which is not yet 
well established. One should keep in mind that different 
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values of the coupling constant for the a- and 7r-bands 
are expected from the calculations^ and observed in the 
dHvA experiment. 7 

With a lack of large single crystals of MgB2, espe- 
cially along the c-axis, the optical measurements were 
done on polycrystalline samples ^iSiSfliSi disorientedS2iS& 
or ab-orientediSiifiii£ films as well as ab-oriented crys- 
tal mosaics li 7 ^ We believe that a large spread of the 
published optical results can be explained by (i) differ- 
ent purity levels of the samples used, (ii) experimental 
difficulties to extract the anisotropic complex conductiv- 
ity from the standard measurements on such objects and 
(iii) a fast rate of the surface degradation in airiSSiS 

In this work, we do optical experiments on individual 
freshly polished single crystals of different purity having 
no contact with the air. The real and imaginary parts of 
the in-plane and out-of-plane optical conductivities are 
derived from the reflectivity and ellipsometry measure- 
ments on the ab- and, for the first time, on the ac-oriented 
crystals. The measurements show a large (« 6.3 eV) and 
almost isotropic plasma frequency. We find different po- 
sitions of the reflectivity plasma edges for the two light 
polarizations. However, such anisotropy is due to extra 
screening of charge carriers by an interband transition at 
2.6 eV, optically active along the ab-plane, rather than 
the anisotropy of the bare plasma frequency itself. The 
anisotropic spectra are in satisfactory agreement with the 
first-principles calculations of the band structure and the 
electron-phonon interaction. We discuss the possible cor- 
rections to the electronic structure which follow from the 
optical results. 

Finally, the presented data tell us what is the 'true' 
color (or rather, colors) of MgB2. 



II. SAMPLES AND PREPARATION 

High-quality single crystals of MgB2 have been grown 
using a cubic anvil technique via the peritectic decompo- 
sition reaction of the MgNBg and Mg at temperature up 
to 2000°C under a pressure of 30 - 35 kbar. The details 
of the crystal growth and extensive characterization are 
given elsewhere*^ 7 . The measurements were done on two 
samples, referenced below as SI and S2. The samples 
were grown under slightly different conditions. For the 
sample SI the purity of Mg was 99.8%, while a 4N-pure 
(99.98%) magnesium was used in the second case. The 
maximum temperature of crystal growth was higher for 
the sample S2 by about 80°C. The time at highest tem- 
perature was 10 minutes for the first sample and 30 min- 
utes for the second one. The T c of both samples is close 
to 38 K (see Fig^), although sample S2 shows somewhat 
narrower transition and a larger field-cooled diamagnetic 
signal compared to the one of SI. 

The dimensions (a x b x c) of the as-grown crystals were 
0.7 x 0.5 x 0.27 mm 3 for SI and 0.6 x 0.5 x 0.18 mm 3 
for S2. We selected a thicker sample SI to prepare the 
(ac) optical face, while sample S2 was used to measure 
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FIG. 1: Magnetization measurements of the superconduct- 
ing transition of the two single crystals of MgB2 used in this 
paper. 



on the (ab) face. The faces were dry-polished using a 0.1 
/im diamond abrasive, since the as-grown surfaces were 
not suitable for quantitative optical examination. 

As it was noticed before^^ the exposed surface of 
MgB2 deteriorates quickly, largely due to the air mois- 
ture. In order to avoid the contamination, the samples 
were kept in a flow of dry nitrogen during and after the 
polishing until the end of the measurements. With this 
precaution, the optical characteristics did not change no- 
ticeably during the experiment. On the other hand, 
switching off the flow immediately affected the optical 
signal as it is shown in the Appendix. 



III. OPTICAL EXPERIMENT AND RESULTS 

Optical properties of MgB2 at room temperature were 
obtained using spectroscopic ellipsometry at 0.75 - 3.7 eV 
and the reflectivity measurements in the infrared range 
from 0.1 to 0.85 eV. In all experiments samples were 
mounted on the top of a sharp optically black cone and 
laser aligned. 

The high-frequency spectra were collected using the 
Woollam VASE32 ellipsometer while the sample was 
kept in a flow of dry nitrogen. The in-plane e a (c<j) = 
ei, a (uj) + Airi(Ji, a {ijj)/oj and the c-axis e c (u>) — ei jC (w) + 
Amtj\ tC {uj)/u) components of the complex dielectric ten- 
sor were both extracted directly from the measurements 
on the ac-surface of the sample SI, using two orthogonal 
crystal orientations and three angles of incidence. For 
the sample S2, only the in-plane optical functions ei, a (w) 
and ai : a(to) were derived from the measurement on the 
ab-plane, using the c-axis data from the sample SI to 
correct for the admixture of the out-of-plane component. 
The details of the recovery of e a {ui) and e c (u>) from the 
ellipsometric output are given in the Appendix. 

The reflectivity R(lu) was measured with a polarizer 
in the range 0.1 - 0.85 eV at a near-normal angle of in- 
cidence using a Fourier transform spectrometer Bruker 
66v. The sample was inside a vacuum chamber of a cryo- 
stat. A gold layer was sputtered in- situ on the crystal 
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FIG. 2: Optical anisotropic spectra of the MgE>2 at 300 K de- 
rived from the ellipsometry and reflectivity measurements on 
the sample SI as described in the text: the normal-incidence 
reflectivity ii(w)(a), optical conductivity oi(u), the dielectric 
function ei(u>) (b), the effective number of carriers N e g (ui) 
and the effective plasma frequency Qp tB g(uj) (c). The a-axis 
and the c-axis spectra are shown by the blue and red colors 
respectively. 



surface to get a reference signal. We reconstructed the 
full reflectivity spectrum using the ellipsometrically de- 
termined dielectric functions in the optical range: 



1 - y/ejjvj 



l+y/e v {Lj) 



, = a,c). 



(1) 



The spectra from two regions were combined in order 
to obtain the complex dielectric function in the whole 
range. While e(ui) was measured directly above 0.75 
eV, at low frequencies we applied a variational Kramers- 
Kronig (KK) constrained analysis of the data. In this 
method, described in Refl28l one finds a KK-consistent 
dielectric function, which gives the best detailed match 
to the reflectivity R(w) at low frequencies and both ei(oj) 



and e2 (w) at high frequencies simultaneously. This pro- 
cedure reduces to a large extent the uncertainty due to 
extrapolations as compared to the usual KK transform. 
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FIG. 3: Comparison of the in-plane optical spectra measured 
on different samples. The same types of spectra are presented 
as in Fig|2] The spectra of samples SI and S2 are shown by the 
solid and dash lines, respectively. Sample S2 has presumably 
a smaller concentration of impurities than sample SI. 

Figs|2t and [2b show the anisotropic reflectivity R(u>), 
optical conductivity o\ (u>) and the dielectric function 
ei(a->), measured on sample SI. Qualitatively, both in- 
plane and c-axis spectra exhibit a similar metallic be- 
havior, characterized by a reflectivity plasma edge, a 
Drude peak in oi(u>) and negative ei(w). However, one 
can see a strong anisotropy. The in-plane reflectivity 
shows a broad plasma edge at about 2 eV. The dielec- 
tric function ei )0 (o;) almost reaches zero at the same fre- 
quency, but then it varies non-monotonically while re- 
maining negative at least up to 3.7 cV. In contrast, the 
c-axis plasma edge is significantly sharper and is higher 
in energy by about 0.5-0.6 eV. Correspondingly, ei iC (a->) 
behaves monotonically and crosses zero at lu* c « 2.6 
eV. These results agree with the previous findings of Fu- 
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damoto and Lee^ 3 . who obtained a plasma edge at about 
16000 cm -1 (2 eV) on a mosaic of ab-oriented crystals. 
Although no direct reflectivity measurements for E \\ c 
were reported so far, the same authors suggested^ 3 , that 
the additional step-like structure at 22000 cm -1 (2.7 eV) 
in the reflectivity of polycrystals samples comes from the 
c-axis plasma edge. The present study fully confirms this 
assignment. 

The in-plane optical conductivity <Tx, a (w) shows an in- 
tense interband peak at ~ 2.6 eV, which is totaly ab- 
sent in the <Ji iC (w) (Fig|2b). Its origin will be discussed 
below. Since this peak is close to the screened plasma 
frequency, it broadens the plasma edge and shifts it to 
lower frequencies by providing an additional screening 
of the charge carriers. It is also responsible for a non- 
monotonic behavior of e\ ta {uj) at higher frequencies and 
even for a structure in i? Q (w) above 3 eV, which resembles 
a 'second' plasma edge. 

Fig. |2t depicts the partial sum rule (effective number 
of carriers) function: 



N e g (to) 



2mV r 



<7i (w')d(jj' , 



(2) 



where m is the free electron mass, V c = 28.9 A 3 is the 
unit cell volume, e is the electron charge, and the corre- 
sponding effective plasma frequency: 



^p,cff(w) 



ai(aj')duj' 



1/2 



(3) 



Since the integration of (J\(lo) has to start from zero fre- 
quency, while the optical data are taken down to 100 
meV, one has to be vigilant about the error bars involved. 
Fortunately, the fact that the dielectric function obtained 
by the method used^ satisfies the KK relations and de- 
scribes both and £2(w) at high frequencies poses 
fairly tight bounds (shown in FigI2t) on the allowed val- 
ues of N e ff . One can see that at the photon energy of 
1.7 eV, which is high enough to comprise most of the 
intraband spectral weight, but below the interband peak 
at 2.6 eV, N e g is about 0.8 for both polarizations, which 
corresponds to the plasma frequency of 6.3 eV. A more 
rigorous estimate which takes into account the broaden- 
ing of the Drude peak, as described below, gives a very 
close value of lu p . This is in contrast with the conclusion 
of Ref 23 about a strong anisotropy of the plasma fre- 
quency, based on the different positions of the reflectivity 
plasma edge for the two polarizations. The anisotropy of 
the reflectivity plasma edge is caused by extra screening 
due to the in-plane interband transition at 2.6 eV and 
not by the strong anisotropy of the unscreened plasma 
frequency. 

It is interesting to compare the in-plane optical proper- 
ties of the two samples (SI and S2), which were prepared 
under different conditions and have presumably slightly 
different impurity levels (see FigEJ ■ O ne can see that the 
two most prominent features of the optical conductivity 
- the Drude peak and the 2.6 eV interband peak - are 



significantly sharper in the sample S2. As a result, the 
reflectivity plasma edge is also sharper and the 'double- 
plasmon' structure in R a (u>) is more pronounced than in 
sample SI. One should keep in mind that the sample S2 
was prepared from a slightly more pure magnesium, and 
it shows a sharper superconducting transition (FigQJ. 
This tells that even a small impurity level (about 0.2% 
in this case) affects significantly optical and transport 
properties. 

The close values of e\, a {uj) for the two samples sug- 
gest that their in- plane plasma frequencies arc similar, 
although the N g is slightly higher for sample SI, which 
could be related to a stronger broadening of the inter- 
band peaks. 



IV. DISCUSSION 

A. Extended Drude analysis 

The extended Drude model has been commonly used 
to analyze interactions in electronic systems. In this for- 
malism, the frequency dependent scattering rate 1/t(cj) 
and effective mass m*(uj)/m are derived from the mea- 
sured complex dielectric function: 
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The only parameters of this conversion are the total 
Drude plasma frequency u p and the high-frequency di- 
electric constant e^, due to all contributions other than 
the conductivity electrons. Both 1/t(u>) and m*(uj)/m 
have direct microscopic interpretation in the context of 
the electron-phonon interaction. One should keep in 
mind, however, that the model assumes that only one 
type of carriers contributes to the Drude response. In the 
case of MgB2, which has two distinct systems of bands 
(cr and 7r), the scattering rate and the effective mass ob- 
tained by Eqs.QJ and (J5J should be regarded as an aver- 
aged value of contributions from each band. Fig0] (sym- 
bols) shows 1/t(w) and m*{u>)/m for the in-plane and 
the c-axis directions, measured on the sample SI. We 
took the values lu pm — 6.28 eV, lu PiC = 6.35 eV, f mi „ = 
11.9, loo^c — 4.77, which are suggested by the data fit- 
ting in Section II V CI (see Table . The value of eoo here 
is given by the sum of in the Table [I] and the oscil- 
lator strengths Si of the Lorentz oscillators. The solid 
curves in Fig0]were calculated using the results of the fit 
described in Section TlV CI 

One can see that the scattering rate and the mass 
renormalization are larger for the in-plane direction. It is 
worth mentioning that the tr-bands with a cylinder-like 
Fermi surface must have a small electromagnetic response 
along the c-axis, while the 7r-bands are expected to have 
comparable contributions in both directions. The present 
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FIG. 4: Extended Drude analysis of the optical conductivity 
of MgB2 (sample SI) at 300 K along the in-plane and the 
c-axis directions. The symbols are the data, the solid curves 
show the two-component fit as described in the Section ft V CI 



result thus suggests that the electron-phonon coupling 
is stronger in the cr-bands. This is in agreement with 
the first-principle calculations of Kong et al., 3 who found 
that the total strength of the electron-phonon interaction 
is dominated by the coupling of the cr-holes to the bond- 
stretching optical phonons. A quantitative analysis of 
the electron-phonon interaction must take into account 
the multi-band electronic structure. 



B. Comparison with ab initio calculations 

A deeper insight can be attained by the comparison 
of the optical data with the existing first-principle cal- 
culations of the band structure and electron-phonon in- 
teraction. One can compute the dielectric function, as- 
suming that the intraband optical conductivity is formed 
by the additive contributions of carriers in the a- and 
the 7r-bands. The two bands are characterized by dif- 
ferent electron-phonon (Eliashberg) transport functions 
a 2 r F(u>) and impurity scattering rates jimp- Intra- 
band optical response can be calculated using finite- 
temperature memory-function formalism for electron- 
boson interaction^ The total conductivity is a sum of 
the intraband, interband (e IB ) and core-electron (e^) re- 
sponses: 



UJ 



p,vf3 



j3—{T,7T 



-lu[lu + iMp(w,T)] 



+ 4 M + e.oo (6) 
(y = a, c). 



The contribution of each of the two bands is determined 
by its respective anisotropic plasma frequency (w Pja , lo p . c ) 



and the memory function^: 

M p (uj, T) = l0imp -2i J dn a 2 tr F p (n)K ( 

where 
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2nT ' 2ttT 
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K(x,y) = - + [tp(l — ix + iy) — ip(l + iy)] + 

y x 



y + x 



[-0(1 -ix- iy) - -0(1 - iy)], 



and ip(x) is a digamma function. 

All the ingredients to compute e(w), except the impu- 
rity scattering rates ja-imp and jnimp, which depend on 
the level and the nature of the impurities (substitutions, 
vacancies, dislocations etc.), are provided by the ab initio 
LDA calculations. 

The spectra of the interband optical conductiv- 
ity <j[ B (u)) obtained by different groups using the 
LMX0 ,30^31 f u ll-potential LAPW— and time-dependent 
DFT— methods, although showing some differences, are 
close to each other. In particular, we took the inter- 
band conductivity spectra, presented in Refl2(il (Fig. 5) 
on the base of the full-potential LAPW method^ and 
computed the corresponding dielectric function e{ B (u>) 
by the KK transformation. The plasma frequencies re- 
ported by different groups^2*2iiS are also close to each 



4.14 eV, 



= 5.89 



other. We used the values w Pi0<T 
eV, uj p , ca = 0.68 eV and lo p .^ = 6.85 eVZLZi 

The ab initio calculated electron-phonon interaction 
functions were presented in Ref3. We used effective 
interaction functions for the two bands a 2 r F a (uj) = 

y 2 r F n ^(uj) + 



a 2 r F aa (uj) + a 2 r F a7r (w) and a&.JVM 



a 2 r F na (u) from the Refs 



The electron-phonon 



coupling constant \ tr — 2 dO,a 2 r F(Q,)/Cl is calculated 
to be about 1.1 for the cr-bands and 0.55 for the 7r-bands. 

The calculated optical spectra are presented in Fig|5] 
Here the scattering rates j a imp — 12.4 meV and ^^imp = 
85.6 meV were taken; this choice is substantiated by the 
data fitting in the next Section. Comparing Figs |21 and [51 
one can notice a remarkable overall agreement and even 
a reasonable quantitative match. 

We begin with the interband transitions. Both theo- 
retical and experimental interband optical conductivities 
show a very strong anisotropy. Notably, the theory pre- 
dicts no sizeable optical intensity along the c-axis of the 
interband transitions below 4 eV, whereas there are two 
low-lying peaks at 0.35 eV and 2.4 eV for polarization 
parallel to the boron planes (see the blue dashed curve 
in FiglSJa). The 2.4 eV peak is due to a transition from 
the cr-band to the 7r-band close to the M-point, where 
a van Hove singularity strongly enhances the density of 
states— The peak at 0.35 eV corresponds to the transi- 
tion between two close cr-bands. The a carriers are local- 
ized within the boron planes, which explains why these 
two transitions are optically very weak along the c-axis. 

Obviously, a strong peak in the experimental spectrum 
along the a-axis at 2.6 eV matches the theoretical peak 
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FIG. 5: Anisotropic optical spectra of MgB2 at 300 K, 
calculated using the results of the published first-principle 
studies-*— as described in the text. The same types of spectra 
are shown as in Fig|5] Additionally, the interband contribu- 
tion along the a-axis is presented in (b) as a dashed blue line. 
The interband conductivity along the c-axis is negligible in 
the shown energy range. 



at 2.4 eV. The discrepancy in the peak position sug- 
gests that the separation between cr-band and 7r-band 
is about 0.2 eV bigger than predicted by the theory. It 
is remarkable^ that the same shift brings the results 
of the dHvA experiments' 1 ' in MgB2 close to theoretical 
predictions 

According to the calculation, the a — > a transition at 
0.35 eV should manifest itself as a noticeable dip in the 
reflectivity spectrum (see inset in Fig|3Ji). Such a dip is 
not observed in the experimental spectra (FigI2K) ■ This 
peak may thus be shifted to even lower frequencies or 
heavily overdamped. The first possibility is consistent 
with a shift of the cr-bands with respect to the Fermi en- 
ergy compared to the calculation, as proposed in Reflssl 
The broadening can be caused by the interaction with the 
conduction electrons, since the energy of this transition 



is within the width of the Drude peak. 

A very intense sharp peak in the out-of-plane conduc- 
tivity a\ tC {uj) is expected around 5 eV i 2 ?^?'^? This excita- 
tion was closely studied in Refl33l by the time-dependent 
DFT method. In the band picture, it comes from a tran- 
sition between almost parallel bands. From the real-space 
point of view, this mode involves charge fluctuations be- 
tween B and Mg sheets, dynamically screened by the in- 
traband transitions. 33 As a consequence, a sharp plasma 
mode at 2.5 eV should emerge. Although 5 eV is be- 
yond our experimental range, there is a sizeable increase 
of f7i. c (w) above 3 eV, which can be a low-energy tail 
of this mode. The c-axis plasmon is at 2.6 eV, which is 
only slightly higher than the calculated value (~ 2.5 eV). 
This is in agreement with the inelastic X-ray scattering 
experiment^ 

The shape and the width of the Drude peak agree well 
with the experiment, which suggests that phonons and 
impurities are the main factors of electron scattering. 
The integrated spectral weight N e g (uj) grows faster ac- 
cording to the calculations, especially for the in plane 
response. This indicates that the theory overestimates 
the value of the plasma frequency. 



C. Two-band fitting of the spectra 

The overall good agreement between the experimental 
and calculated data allows us to use the chosen model to 
fit the spectra, treating the plasma frequencies and scat- 
tering rates as adjustable parameters. We could achieve 
a satisfactory least-square fit of reflectivity R(uj) in the 
mid-infrared and ei(w) and <j\(uj) at higher frequencies 
using formulas © and JJJ) as it is shown in Fig El The 
corresponding parameter values are given in Tabled 

It turns out that leaving all four plasma frequencies ad- 
justable makes the fitting procedure under-determined. 
Therefore, assuming the 2D nature of the cr-band, we 
fixed the plasma frequencies of the cr-band to the values 
given by ab initio calculations, and left only the total in- 



plane and c-axis plasma frequencies cop 1 tot 
adjustable. Another assumption was that the ab initio 
calculations^iSiSi correctly describe the electron-phonon 
functions oq r F atV (w). 

We modelled the interband conductivities by Lorentz 
oscillators (v — a, c): 



JB 



(-) = £■ 



(8) 



with adjustable frequency to V i, oscillator strength S U i and 
width 7^. Keeping in mind the two interband peaks be- 
low 3 eV predicted by the theory, we put two Lorentzians 
to model £q S (w). Only one oscillator term above 3 eV 
was taken for e^ B (w). 

One can see that the bare plasma frequencies w PiCr and 
uip^ are almost the same and equal to 6.3 eV, which 
confirms our previous estimate based on the partial sum 
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FIG. 6: Experimental data (symbols) and the multi-band fit 
(solid lines) of the reflectivity (a) and the dielectric function 
(b) along the a-axis (blue) and the c-axis (red). 



rule (Fig|2t)- This value is much higher compared to 
previous reporto 14 i 19 i 20 i 21 i 22 i 23 of 1.5-2.5 eV. Thus, the 
discrepancy with the theoretical value of 7 eV is likely 
to be much less than it was thought before. However, 
the current mismatch is not negligible since it results in 
about 20-25% deviation of the Drude spectral weight. 
It is worth to mention that the extremal orbit areas in 
the de Haas-van Alphen experiment* on both a and it 
Fermi surfaces are also somewhat smaller than predicted 
by theory It was pointed outjiS that the discrepancy 
can be removed by a shift of the cr-bands downward by 
about 115 meV and the 7r-bands upward by 125 meV. 
This is in perfect agreement with our observation of the 
mismatch of 0.2 eV in the position of the 2.6 eV peak, as 
it is mentioned before. This also qualitatively explains 
the smaller value of the plasma frequency. 

The impurity scattering rate in the 7r-band y-nimp is 
about 85 meV. Since the 7r-electrons have rather strong 
optical spectral weight and a modest electron-phonon in- 
teraction, the absolute values of both R a (oj) and R c (w) 
are quite sensitive to this parameter. In contrast, the 
value of jaimp does not significantly affect the spec- 
tra and cannot be accurately determined from the fit, 
because a large electron-phonon interaction in the a- 
band dominates the total scattering above 100 meV at 
300 K. However, y a imp influences drastically the shape 
of the temperature-dependent in-plane DC resistivity 
p a {T) = l/(Ti_ a (uj — > 0, T) curves, most notably the value 
of p(300K)/p(40K) (RRR), as can be seen from Figd 



TABLE I: Model parameters which give the best match to 
the measured optical spectra of MgB2 at 300 K and the ra- 
tio p a (300K) / p a (40K) from Refl37l as described in the text. 
Parameters, marked by *, were fixed to the values given by 
the band calculations. 



Parameter ab-plane c-axis 





4.14* eV 0.68* eV 




4.72 eV 6.31 eV 




6.28 eV 6.35 eV 


r yaimp 


12.4 meV 


"7 V imp 


85.6 meV 


LU 1 


2.60 eV 3.92 eV 


Si 


1.79 1.7 


7i 


1.11 eV 1.51 eV 


Ul2 


8.57 eV 


S 2 


6.81 


72 


100.2 eV 




3.31 3.07 



A measurement*^ on a crystal grown in the same group 
under similar condition, showed RRR ss 5 along the a- 
axis. Therefore we have chosen of Jaimp = 12.4 meV, 
which gives the same resistivity ratio (see Fig[7|). The as- 
sumption "faimp = y-Kimp gives RRR less than two, which 
strongly suggests that y a imp is actually several times 
smaller than ^imp- This agrees with the Raman study^ 
on single crystals, where a relation jTrimp/jaimp = 6 — 9 
was found, although the absolute values of the scattering 
rates are smaller than in our case. 




100 200 300 
T(K) 



100 200 300 
T(K) 



FIG. 7: Temperature-dependent resistivity curves along the 
a-axis (blue) and c-axis (red), calculated using the results of 
the fit of the optical data for 7^i mp = 85.6 meV and different 
values of the scattering rate in the er-band: 7 CT i mp = 12.4 meV 
< 7ir™ P (left) and y aimv = 7^mp (right). 

A disparity between the impurity scattering rates in 
the two bands and a small interband a <-> ir scattering 
were proposed in Refl2ll to explain the surprisingly small 
dependence of T c on the impurity level, not expected for 
a two-band superconductor. Microscopically, this can be 
explained by the fact that the electronic wavefunctions 
in cr-bands are confinedto the boron planes and not ef- 
ficiently scattered by the magnesium vacancies and sub- 
stitutions. On the other hand, the same defects strongly 
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FIG. 8: Two images of the same sample (ac-plane) made in 
a polarized light: E || ab (left) and E \\ c (right). The sample 
looks golden for E \\ c, due to a sharp reflectivity plasma edge 
at around 2.5 eV, while the plasma edge is smeared for E \\ ab 
by a strong a — > it interband transition at 2.6 eV. 

scatter the electron states in the 3D 7r-bands. Our results 
support this conjecture. 

The first Lorentz term for the a-axis at 2.6 eV clearly 
corresponds to the discussed already interband transi- 
tion. However, the second oscillator is extremely broad 
and it cannot be matched with a predicted interband 
peak at 0.35 eV. This probably means this relatively weak 
peak is shifted and/or broadened so much that it cannot 
be identified in the optical spectrum. The second Lorentz 
term may thus describe the background formed by the 
tails of several broadened interband peaks. 

D. The color of the magnesium diboride 

Finally we address a practical question: what is the 
color of MgB2? The characterizations given in the liter- 
ature varied significantly, ranging from the 'golden' and 
'tan' to 'silver-metallic', 'black' and even 'blue'. Lea22, 
has noticed that, depending on the polarization, "...the 
fresh surface of as-grown crystals in the ab-plane changes 
the color from silver to dark gray. In contrast, for the 
ac-plane of the crystals a beautiful change of color from 
golden- yellow to blue can be seen." Fig|8] shows the two 
images of the unpolished ac-plane of the sample SI for 
different polarizations of the incident white light. One 
can see that the sample spectacularly changes color from 
a blueish-silver for E || ab to the yellow for E \\ c. The 
explanation comes naturally from the reflectivity curves 
shown in Fig[5]i. For the c-axis polarization, it is a sharp 
plasma edge at about 2.5 eV, the same as of gold, which 
makes the sample yellow. For the ab-plane polarization, 
the plasma edge is smeared due to the strong interband 
transition at 2.6 eV. As a result, the reflectivity spectrum 
in the visible range (1.8 - 3.1 eV) is relatively flat with 
a maximum (stronger in sample S2) at ~ 2.8 eV which 
makes the color blue-silver. 

Knowing the anisotropic spectra, we can predict the 
reflectivity shape of MgB2 in the polycrystalline form. 
If the light wavelength A is much larger than a typical 
grain size I, one can apply the effective- medium approx- 
imation (EMA) as prescribed in Reflifl. In the short- 



wavelength limit A <C I, it is more adequate to aver- 
age directly the reflectivities along the two axes, as it 
was done in Refill In FigEl both calculated curves 
Rema{u) (for sphere-like and plate-like crystallites) and 
R av {uo) = (2/3)i? a fc(cj) + (1/3)R c (uj) are shown, together 
with the spectrum on a polycrystal from our previous 
study iSi Notably, Rema{oj) and i? a „(w) are close to each 
other and their overall shape in the visible range matches 
well the measured spectrum. The latter shows only a 
weak step-like feature near 2.5 eV, coming from the o 
axis contribution^ which explains why the polycrys- 
talline samples have typically black, or slightly tan, color. 
The absolute value of the experimental reflectivity from 
Ref I20I is by 10-20% lower than the calculated one, which 
is likely caused by an overdamping of the 7r-bands due to 
a much higher impurity level (presumably MgO) of the 
sample used. 

Energy (eV) 
12 3 




0.2 - 



0.0 1 ' ' ' ' • 1 

10000 20000 30000 

Wavenumber (cm 1 ) 

FIG. 9: Reflectivity spectra of polycrystalline MgB2 at 300 
K. The spectra, calculated using the effective-medium ap- 
proximation (spherical or plate-like crystallites) and the di- 
rect reflectivity averaging, are compared with experimental 
spectrum from Refl2bl 



V. CONCLUSIONS 

We investigated anisotropic optical properties of MgB 2 
single crystals of different purity at room temperature in 
the energy range 0.1 - 3.7 eV and compared them with 
the existing first principle calculations of the electronic 
structure and electron-phonon coupling. The analysis of 
the anisotropy allowed us to distinguish properties of the 
a- and 7r-bands. 

The total bare plasma frequencies along the in-plane 
and c-axis directions are almost the same and equal to 6.3 
eV, which is much higher than the previously reported 
values. However, it is still smaller than the theoretical 
value of about 7 eV. 

The shape of the Drude peak is well described by the 
electron scattering on phonons and impurities. The data 
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suggest that the impurity scattering in the 7r-bands is sev- 
eral times larger than in the cr-bands in agreement with a 
proposal of Refill! aimed to explain a surprisingly small 
dependence of the critical temperature on the sample pu- 
rity. The electron-phonon interaction is stronger in the 
cr-bands, in agreement with the theory^ and de Haas- van 
Alphen experiments^ 

The in-plane optical conductivity clearly shows an in- 
tense peak at 2.6 eV, which corresponds to a transition 
from the a- to the 7r-band. It is higher by about 0.2 eV 
than the theoretical value. An interband peak due to the 
transitions between two close cr-bands, expected to be 
at 0.35 eV, was not seen in the experiments which can 
be explained by an overdamping or a shift towards lower 
frequencies relative to the theoretical calculations. 

The width of the Drude peak and the interband a — > tt 
peak at 2.6 eV are very sensitive to the sample purity. 
This can partially explain the spread of the results of 
previous optical studies. 

The two quantitative mismatches between the exper- 
iment and the theory, namely (i) a smaller plasma fre- 
quency, (ii) a higher by 0.2 eV energy of the 2.6 eV 
cr — ► 7r transition, tell us that the cr-bands are proba- 
bly shifted down and 7r-bands are shifted up compared 
to the calculations, so that their relative shift is about 0.2 
eV. Notably, the same conclusion was derived in Refl35l 
from the analysis of the dHvA measurements £ 

The positions of the reflectivity plasma edges for two 
polarizations are very different (~ 2 eV for a-axis and ~ 
2.5 eV for the c-axis). This is caused by an additional 
screening of the intraband carriers by the interband tran- 
sition at 2.6 eV, but not due to different ab-plane and 
c-axis bare plasma frequencies. As a result, the color of 
the sample depends on the polarization of the light: it is 
blueish-silver for E || a and yellow for E |j c. 
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MaNEP. We are grateful to O.V. Dolgov, J. Kortus and 
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Appendix. Determination of e a (uj) and e c (u)) from 
ellipsometry 

The general geometry of the ellipsometric experiment 
in the configuration 'fixed analyzer (at 45°) - sample - 
rotating analyzer' is shown on FigllOb. The axes x, y 
and z are assumed to be along the principal axes of the 
dielectric tensor of the sample. Ellipsometry provides 
two parameters ip and A, related to the ratio of reflection 
coefficients for the p- and s-polarizcd light: 

p = r p /r s = tan^)exp(iA). (9) 

r p and r s are given by the Fresnel formulas: 

_ VT - Sz 1 sin 2 9 - ^fsZcosQ 
Vl-e z sin 9 + ^feZ cos 9 

cos 9 — \/e v — sin 2 9 

r s = V ■ ' 

cos 9 + J e y — sin 2 9 

where e„ = e\ v + itiv (y — x, y, z) are the components 
of the complex dielectric tensor. In the case of optically 
uniaxial MgB-2 (e a = f-b =/= £e) three different orientations 
(xyz) are possible: (aac), (aca) and (caa) (see FigllOb). 




(aac) (aca) (caa) 



FIG. 10: (a) The general configuration of ellipsometric ex- 
periment; (b) three specific sample orientations used in this 
paper. 



On the sample SI we did ellipsometry in the orien- 
tations (caa) and (aca). This yields four independent 
quantities i/w, A aca , ip caa and A caa at every photon 
energy and every chosen angle of incidence. Since each 
of these functions depend on the four values ei. a , £2. a, ci.e 
and £2,cj the latter ones can be obtained by the numerical 
inversion of the four Fresnel equations. This procedure, 
when applied to three different angles of incidence (60°, 
70° and 80°), gave close results. In order to improve the 
accuracy of the output, we determined by the least square 
fitting the values of e±. a , £2,0., ci.e and e2 jC that render the 
best match to the measured ellipsometric parameters at 
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all mentioned angles of incidence simultaneously. The 
experimental and fitting curves for the sample SI for the 
orientations (aca) and (caa) are shown in Figs lllh , and 
ITTb respectively. 
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FIG. 11: Ellipsometric spectra t[> and A for three angles of 
incidence (60°, 70°, and 80°) in the orientations (aca) (a) 
(caa)(b) and (see Fig llOb 1 ). measured on sample SI. Sym- 
bols are the measurement results, the solid line correspond to 
the fit as described in the text. 

The second experiment was done on the sample S2 in 
the orientation (aac) . The corresponding spectra of ipaac 
and A ooc for the same three angles of incidence are shown 
in Fig^](as solid symbols). These data, taken alone, are 
not sufficient to extract both e a and e c independently. 
Since ip aac and A aac are not very sensitive to the value 
of e c except close to the screened plasma frequency (ac- 



cording to Aspnes approximation)^ we extracted e a (u>), 
assuming that e c (w) is the same as the sample SI. One 
should realize that this introduces some uncertainty, be- 
cause the samples SI and S2 might have somewhat dif- 
ferent c-axis dielectric functions. This uncertainty is neg- 
ligible below 2 eV, but somewhat enhanced close to the 
c-axis plasmon at 2.6 eV. The best fit of if> and A is given 
by the solid lines in Fig|P2l 
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20000 
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30000 



FIG. 12: Ellipsometric spectra tp and A for three angles of 
incidence (60°, 70°, and 80°) in the orientation (aac) (see 
FigllOb') . measured on sample S2. Symbols are the measure- 
ment results, the solid line correspond to the fit as described 
in the text. 

Fig ll3l shows how quickly the ellipsometric parameters 
change when a surface is exposed to the air, due to the 
formation of a contamination layer. 




1000 2000 

Time (sec) 

FIG. 13: The effect of the exposure to the air on the optical 
properties of MgB2. The time dependence of ellipsometric 
parameter A at 1.5 eV is shown (sample SI, geometry (caa), 
angle of incidence 70°). Initially the sample was kept in a 
flow of dry nitrogen, which was switched off at the moment 
designated by the arrow. 



11 



1 J. Nagamatsu, N. Nakagawa, T. Muranaka, Y. Zenitani, 
and J. Akimitsu, Nature 410, 63 (2001). 

2 I. I. Mazin and V.P. Antropov, Physica C 385, 49 (2003). 

3 Y. Kong, O.V. Dolgov, O. Jepsen, and O.K. Andersen, 
Phys. Rev. B 64, 020501 (R) (2001). 

4 S. L. Bud'ko, G. Lapertot, C. Petrovic, C.E. Cunningham, 
N. Anderson, and P.C. Canfield, Phys. Rev. Lett. 86, 1877 
(2001). 

5 D.G. Hinks, H. Claus, and J.D. Jorgensen, Nature 411, 
6836 (2001). 

6 H. Uchiyama, K.M. Shen, S. Lee, A. Damascelli, D.H. Lu, 
D.L. Feng, Z.-X. Shen, and S. Tajima, Phys. Rev. Lett. 
88, 157002 (2002). 

7 E.A. Yelland, J.R. Cooper, A. Carrington, N.E. Hussey, 
P.J. Meeson, S. Lee, A. Yamamoto, and S. Tajima, Phys. 
Rev. Lett. 88, 217002 (2002). 

8 Y. Wang, F. Bouquet, I. Sheikin, P. Toulemonde, B. Revaz, 
M. Eisterer, H.W. Weber, J. Hinderer, and A. Junod, J. 
Phys.: Condens. Matter 15, 883 (2003) . 

9 F. Bouquet, Y. Wang, I. Sheikin, T. Plackowski, A. Junod, 
S. Lee, and S. Tajima, Phys. Rev. Lett. 89, 257001 (2002). 

10 A. A. Golubov, J. Kortus, O. V. Dolgov, O. Jepsen, Y. 
Kong, O. K. Andersen, B. J. Gibson, K. Ahn and R. K. 
Kremer J. Phys.: Condens. Matter 14, 1353-1360 (2002). 

11 M. R. Eskildsen, M. Kugler, S. Tanaka, J. Jun, S. M. Kaza- 
kov, J. Karpinski, and O. Fischer, Phys. Rev. Lett. 89, 
187003 (2002). 

12 S. Galambosi, J. A. Soininen, A. Mattila, S. Huotari, S. 
Manninen, Gy. Vanko, N. D. Zhigadlo, J. Karpinski, and 
K. Hamalainen, Phys. Rev. B. 71, 060504(R) (2005). 

13 B. Gorshunov, C.A. Kuntscher, P. Haas, M. Dressel, F.P. 
Mena, A.B. Kuzmenko, D. van der Marel, T. Muranaka, 
and J. Akimitsu, Eur. Phys. J. B 21, 159 (2001). 

14 R.A. Kaindl, M.A. Carnahan, J. Orenstein, D.S. Chemla, 
H.M. Christen, Hong-Ying Zhai, M. Paranthaman, and 
D.H. Lowndes, Phys. Rev. Lett. 88, 027003 (2002). 

15 A. Pimenov, A. Loidl, and S.I. Krasnosvobodtsev, Phys. 
Rev. B 65, 172502 (2002). 

16 J.H. Jung, K.W. Kim, H.J. Lee, M.W. Kim, T.W. Noh, 
W.N. Kang, Hyeong-Jin Kim, Eun-Mi Choi, C.U. Jung, 
and Sung-Ik Lee, Phys. Rev. B 65, 052413 (2002). 

17 A. Perucchi, L. Degiorgi, J. Jun, M. Angst, and J. Karpin- 
ski, Phys. Rev. Lett. 89 097001 (2002). 

18 M. Ortolani, D. Di Castro, P. Postorino, I. Pallecchi, M. 
Monni, M. Putti, and P. Dore, Phys. Rev. B 71, 172508 
(2005). 

19 J.J. Tu, G.L. Carr, V. Perebeinos, C.C. Homes, M. Stron- 
gin, P.B. Allen, W.N. Kang, Eun-Mi Choi, Hyeong-Jin 



Kim, and Sung-Ik Lee, Phys. Rev. Lett. 87, 277001 (2001). 

20 A.B. Kuzmenko, F.P. Mena, H.J. A. Molegraaf, D. van der 
Marel, B. Gorshunov, M. Dressel, I.I. Mazin, J. Kortus, 
O.V. Dolgov, T. Muranaka, and J. Akimitsu, Solid State 
Commun. 121, 479 (2002). 

21 I.I. Mazin, O.K. Andersen, O. Jepsen, O.V. Dolgov, J. 
Kortus, A. A. Golubov, A.B. Kuzmenko, and D. van der 
Marel, Phys. Rev. Lett. 89, 107002 (2002). 

22 Mi-Ock Mun, Young Jin Kim, Yu Park, Jae Hoon Kim, 
S.H. Moon, H.N. Lee, H.G. Kim, and B. Oh, J. of Super- 
cond. 15, 475 (2002). 

23 Y. Fudamoto and S. Lee, Phys. Rev. B 68, 184514 (2003). 

24 F. Marsiglio, Phys. Rev. Lett. 87, 247001 (2001). 

25 E.G. Maksimov, J. Kortus, O.V. Dolgov, and I. I. Mazin, 
Phys. Rev. Lett. 89, 129703 (2002). 

26 D. Chvostova, V. Zelezny, L. Pajasova, A. Tarasenko, A. 
Plecenik, P. Kus, and L. Satrapinsky, Thin Solid Films 
455, 213 (2004). 

27 J. Karpinski, M. Angst, J. Jun, S. M. Kazakov, R. Puz- 
niak, A. Wisniewski, J. Roos, H. Keller, A. Perucchi, L. 
Degiorgi, M. R. Eskildsen, P. Bordet, L. Vinnikov and A. 
Mironov, Supercond. Sci. Technol. 16, 221 (2002). 

28 A.B.Kuzmenko, Rev. Sci. Instrum. 76, 083108 (2005). 

29 O.V.Dolgov and S.V. Shulga, J. Supercond. 8, 511 (1995). 

30 V.P. Antropov, K.D. Belashchenko, M. van Schilfgaarde, 
and S.N. Rashkeev, cond-mat/0107123 

31 P. Ravindran, P. Vajeeston, R. Vidya, A. Kjekshus, and 
H. Fjellvag, Phys. Rev. B. 64 224509 (2001). 

32 J. Kortus, I.I. Mazin, K.D. Belashchenko, V.P. Antropov, 
and L.L. Boyer, Phys. Rev. Lett. 86,4656 (2001). 

33 W. Ku, W.E. Pickett, R.T. Scalettar, and A.G. Eguiluz, 
Phys. Rev. Lett. 88, 057001 (2002). 

34 O.V.Dolgov, Y. Kong, J. Kortus and O. Jepsen, private 
communication. 

35 H. Rosner, J.M. An, W.E. Pickett, and S.-L. Drechsler, 
Phys. Rev. B 66, 024521 (2002). 

38 I.I. Mazin and J. Kortus, Phys. Rev. B, 65 180510(R) 
(2002). 

37 A.V. Sologubenko, J. Jun, S.M. Kazakov, J. Karpinski, 
and H.R. Ott, Phys. Rev. B 66, 014504 (2002). 

38 J.W. Quilty, S. Lee, S. Tajima, and A. Yamanaka, Phys. 
Rev. Lett. 90, 207006 (2003). 

39 S. Lee, Physica C. 385 31 (2003). 

40 RE. Sulewski, T.W. Noh, J.T. McWhirter, and A.J. Siev- 
ers, Phys. Rev. B 36, R5735 (1987). 

41 D.E. Aspnes, J. Opt. Soc. Am. 70, 1275 (1980). 



